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A dielectric drop suspended in conducting liquid and subjected to an uniform electric field deforms 
into an ellipsoid whose major axis is either perpendicular or tilted (due to Quincke rotation effect) 
relative to the applied field. We experimentally study the effect of surface-adsorbed colloidal particles 
on these classic electrohydrodynamic phenomena. We observe that at high surface coverage (> 90%), 
the electrohydro dynamic flow is suppressed, oblate drop deformation is enhanced, and the threshold 
for tilt is decreased compared to the particle-free drop. The deformation data are well explained 
by a capsule model, which assumes that the particle monolayer acts as an elastic interface. The 
reduction of the threshold field for rotation is likely related to drop asphericity. 

PACS numbers: 


I. INTRODUCTION 


Electric fields provide a versatile means to control 
small-scale fluid and particle dynamics, e.g., electrohy¬ 
drodynamic instabilities for pattern formation in thin 
polymer films m or particle suspensions j4], electro¬ 
hydrodynamic atomization to produce micro- and nano¬ 
particles mi, drop and vesicle manipulation m hb> 
and colloidal assembly HHH]. Fluid interfaces pro¬ 
vide additional functionality opening new routes for the 
bottom-up fabrication of novel structurally complex ma¬ 
terials p~01 [22] . Recent works [231426] find that micro¬ 
particles constrained on a drop surface can form vari¬ 
ous structures in the presence of applied uniform electric 
field. The underlying mechanisms are still under investi¬ 
gation but a major driving force in this system is the flow 
created by the electric shear stresses due to accumulation 
of charges at the interface. 

A particle-free drop placed in an electric field polar¬ 
izes because of the mismatch of the bulk fluids electrical 
conductivity, cr, and dielectric constant, 6 
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Upon application of an electric field, mobile charges car¬ 
ried by conduction accumulate at the boundary (even 
though the net charge on the interface remains zero), see 
Figure [l] For a sphere placed in a uniform electric field 
with magnitude Eq, the induced surface free charge dis¬ 
tribution increases as [27, EB, 
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FIG. 1: Drop polarization and electrohydrodynamic flow stream¬ 
lines for RS < 1 in applied uniform DC electric field, (a) Weak 
fields with Eq < Eq , where Eq is given by Eq. [7] induce pure 
straining flow and axisymmetric oblate deformation, (b) In strong 
fields, Eq > Eq constant torque is induced by the misaligned 
dipole, the flow acquires a rotational component, and the drop 
is tilted with respect to the applied field direction. 


where 0 is the angle with the applied field direction and 
t mw is the Maxwell-Wagner polarization time. The sub¬ 
scripts “in” and “ex” denote the values for drop and sus¬ 
pending medium, respectively. Thus, the polarity of the 
induced charge is determined by the product of R and 
S , which compares the conduction response of the me¬ 
dia [29, [30] . If RS < 1, the conduction in the drop is 
slower than the suspending liquid. As a result, the inter¬ 
face charge distribution is dominated by charges brought 
from the suspending medium and the drop dipole mo¬ 
ment is oriented opposite to the applied electric field di¬ 
rection (note that in the opposite case, RS > 1, the par¬ 
ticle dipole is reversed and it is aligned with the electric 
field.) 

The electric field acting on the induced surface charges 
creates a tangential stress, r r e = QEq. In the case of a 
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simple fluid-fluid interface such as oil-water this shear 
stress can only be balanced by viscous stresses due to 
fluid flow. In weak fields, for which the electric stresses 
are smaller than the capillary stress due to surface 
tension 7 /a so that drop shape remains nearly spherical 
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where a is the drop radius and 7 is the interfacial ten¬ 
sion, at steady state the fluid undergoes axisymmetric 
straining flow about the drop m , see Figure [l] (a), with 
surface velocity 
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where A = /ii n /p ex is the viscosity ratio. The charge dis¬ 
tribution corresponding to RS < 1, illustrated in Figure 
[I] (a), results in surface fluid motion from the pole to the 
equator and the equator is a stagnation line. Drop defor¬ 
mation resulting from the action of the electric stresses 
is either oblate or prolate ellipsoid EH 
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where d\\ and d± are the spheroid axes parallel and per¬ 
pendicular to the direction of the applied electric field. 

Recently it was observed that a nonaxisymmetric rota¬ 
tional flow may appear for drops with RS < 1, see Figure 
[T] (b). The drop can assume steady tilted orientation with 
respect to the electric field direction [32, 33] or oscillate 
IS, 35 . These behaviors have been linked to the Quincke 
rotation phenomenon [36H39] , which is an instability aris¬ 
ing from the unfavorable orientation of the induced dipole 
(for RS < 1 the dipole direction is opposite to the applied 
electric field and “wants” to flip). The complete flip of 
the dipole is prevented by charge supply from the bulk: 
the induced surface charge (dipole) rotates with the par¬ 
ticle, but at the same time the suspending fluid recharges 
the interface. The balance between charge convection by 
rotation and supply by conduction from the bulk results 
in a steady misaligned torque and continuous spinning 
in the case of a rigid sphere, the so called Quincke rota¬ 
tion. The rotation rate uo and the steady oblique dipole 
orientation, characterized by the angle a (illustrated in 
Figure [l](b)) are 
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where Eq is the threshold electric field above which the 
rotation occurs and t mw is the Maxwell-Wagner polar¬ 
ization time 
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R and S are the conductivity and permittivity ratios de¬ 
fined by Eq. [l] In the case of drops, the electric torque 
generates rotational fluid motion [32, 3T, in addition to 
the already present straining flow. The resulting linear 
flow causes the drop to deform into a general ellipsoid 
whose major axis is misaligned with the applied electric 
field.The steady oblique orientation and deformation of 
the drop were experimentally and theoretically investi¬ 
gated by our group [33, (35, 3D] . 

Given the complex behavior of a drop with a simple 
fluid interface, a question naturally arises: how does sur¬ 
face modification changes drop electrohydrodynamics? 
We studied experimentally a microparticle-coated drop 
with RS < 1 in the Taylor regime (below the thresh¬ 
old for Quincke rotation) and at low to moderate sur¬ 
face coverages (below 50%) [24] . Particles initially ran¬ 
domly distributed at the surface formed a “belt” around 
the equator as the drop deforms into an oblate shape. 
This is an expected consequence of the straining elec¬ 
trohydrodynamic flow. In stronger DC fields, the dy¬ 
namics of these “armored” drops becomes more com¬ 
plex [23] |24| and depends on the particles characteris¬ 
tics. Belts formed by low polarizability particles break 
into a sequence of counter-rotating vortices of particles. 
When dipole-dipole attraction becomes strong the par¬ 
ticle chain and the drop experiences a prolate deforma¬ 
tion and tip-streaming occur with ejection of particles. 
For non-spherical conductive particles, we have observed 
drop “kayaking”: its major axis precessing around the 
applied field direction. 

In this work we investigate drops at high particle- 
coverage. The particle-covered drop is modeled as a cap¬ 
sule and the experimentally measured drop deformation, 
D, as a function of field strength, Ca, is used to find 
the shear elasticity of particle monolayers. We also in¬ 
vestigate the effect os particle coverage on the onset for 
Quincke rotation, and we find that high surface coverage 
of particles drastically reduces the critical fields strength 
of rotation. 


II. EXPERIMENTAL SECTION 
A. Materials 

The suspending fluid is castor oil (Alfa Aesar) with 
viscosity p ex = 0.69 Pa.s, dielectric constant £ ex = 4. 6 ^ 0 , 
where So is the vacuum permittivity, conductivity cr ex = 
3.8 x 10 - 11 S/m, and density p ex = 962 kg/m~ 3 . The 
drop fluid is silicon oil with /q n = 0.05 Pa.s, S[ n = 2. 8 ^ 0 , 
<7 in = 2.4 x 10 - 12 S/m, p in = 963.5 kg/m~ 3 (UCT). The 
surface tension between the castor oil and the Silicone 
oil is 4.5 mN/m |33j. The permittivity and the conduc¬ 
tivity were measured using a dielectric constant meter 
and a conductivity meter from Scientifica. The viscosity 
and density are specified by the company-provider. The 
characteristics of the particles are summarized in Table 
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Shape 

Type 

density p v (kg/m' 

3 ) radius r (pm) Conductivity 

Supplier 

irregular 

Aluminum (Al) 

2600 

1.5, 12 

++ 

Atlantic equipment 

sphere 

glass (G) 

2200 

3.5, 8.5 

+ 

Corpuscular/Cospheric 

sphere 

polyethylene (Pe) 

1000 

50 

- 

Cospheric 


TABLE I: Particles characteristics. The symbol + indicates that 
the particles are slightly more conducting that the fluids, ++ highly 
more conducting, and — less conducting. 


B. Methods 



FIG. 2: (a) Drop observation for the experimental measurements, 
(b) D(Ca) for a particle-free drop is used to validate the set-up. 
The symbols are experimental data and the line is the Taylor’s 
theory Eq. [ 5 ] (c) Definition of drop deformation and orientation 
characteristics: the spheroid axes and d± and the tilt (3. 


D = (d|| — d±_)/(d\ \ +d±) and the angle (jr/2 — f3) between 
the drop major axis and the applied field direction, see 
Figure [2] for definitions. The precision of the measure¬ 
ments is set by the resolution of the pictures; for a drop 
of radius a = 2 mm it is around 0.0055 mm/px leading to 
an absolute error in the deformation (D) of 0.0045, which 
decreases for smaller drops. The surface concentration of 
particles, ip, is defined as the percentage of the drop sur¬ 
face covered by particles (including the space between 
particles) once the particles are brought together. We 
measure p at very low electrical fields, Eq < 60 kV/m, to 
minimize particle compaction. The experimental proce¬ 
dure and the material properties are validated by com¬ 
paring drop deformation with the Taylor model for small 
deformat ions [31 j given by Eq. [ 5 ]. 


The experimental set-up has the same design as the 
one used in m- A uniform DC electric field is created in 
a parallel-plate chamber made of two 7.6cm by 10.5 cm 
brass electrodes separated by a gap of 4.5cm. Fields up 
to 16 kV/cm are generated using a voltage amplifier con¬ 
nected to a DC Power supply. The chamber is filled with 
castor oil. A drop of silicon oil with suspended micro¬ 
particles is injected in the middle of the chamber. The 
particles are driven to the interface by the application 
of an electric pulse the duration and strength of which 
depend on the type of micro-particles; after the field is 
turned off the particles remain trapped at the interface as 
the thermal energy is negligible compared to the energy 
required to remove the particle from the interface.The 
drop is then manually moved around in order to ran¬ 
domize the particle distribution on the drop surface. In 
the actual experiment, electric field is applied and drop 
behavior is recorded for about 3 minutes. On this time 
scale drop sedimentation is negligible. After each record¬ 
ing, the electrical field is turned off and the drop is moved 
back to its initial position. This action removes any parti¬ 
cle structures. The experiment is repeated for a different 
electric field strength. 

The drop is observed from a direction perpendicular 
to the field. The top view of the drop is recorded every 
0.2s, which is a compromise between memory limitation 
of the software and the length of the recording. The 
images are post-processed with ImageJ and Matlab to 
extract the drop diameter d, the deformation parameter 


III. RESULTS 

Figure [3] illustrates the typical behaviors of drops with 
p rsj 90% surface coverage. In weak fields, the drop al¬ 
ways deforms into an oblate spheroid and the particles 
accumulate at the equator forming a “belt”. At high cov¬ 
erages, the belt is very wide and only a small region near 
the poles remains particle free. In stronger fields, the 
oblate deformation increases and the drops may adopt 
peculiar “drum-like” shapes. In even stronger fields, the 
drops may start rotating or implode. 



FIG. 3: Typical drop behaviors at high surface coverage of glass 
spheres ( with radius r = 8.5/im). The scale bar corresponds to 
1 mm. 
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A. Weak DC electric fields: increased oblate 
deformation 



FIG. 4: Drop deformation D as a function of the dimensionless 
field strength Ca. The symbol X corresponds to a “clean” (particle- 
free) drop with radius a = 2.2mm, the symbols o, □, o and A denote 
drops covered with A1 (r = 12 fim) with ip = 65%, 88%, 93% and 
96% and a = 1.9mm, 1.8mm, 1.7mm and 2.2mm, respectively. The 
symbol ★ corresponds to a drop covered with A1 (r = 1.5 fim) with 
(p = 80% and a = 2.3mm, the symbol <1 corresponds to a drop 
covered with Pe (cp = 78%, a = 2.2mm), the symbols +, and light 
X denote drops covered with G (r = 8.5 fim) with (p = 100% and 
93% and a = 1.2mm, and 2mm, the symbol > corresponds to a drop 
covered with G (r = 3.5 /am) with cp — 100% and a = 1.7mm. The 
solid line corresponds to the Taylor prediction Eq. [5] for viscosity 
ratio A = 0.07. The inset shows the deformation of a particle-free 
drop in a wider range of Ca. Filled symbols denote drops in the 
electrorotation (tilted) regime. For clarity, error bars are shown 
only for some drops as they are similar for all the measurements. 


Figure [4] shows that particle-covered drops undergo 
larger deformation compared to uncoated (“clean”) drops 
at the same field strength. Moreover, the magnitude of 
the deformation is insensitive to particle type, size, and 
coverage (above ip = 65%). 



FIG. 5: Streamlines of the flow around a drop, a) Theoretical 
flow around a clean drop, b) Trajectories measured around a glass 
covered drop 


Once the drop reaches a steady oblate deformation, 
the particles at the surface do not move, which suggests 
the absence of electrohydrodynamic flow. To test for 



FIG. 6: (a) - Maximum inclination fimax measured over 4 minutes 
as a function of E/Eq for a drop covered with G spheres (r = 
3.5 /am). The symbol + represents the second measurement realized 
for some E/Eq. b- Examples of some temporal evolution of (3. The 
symbol □ corresponds to E/Eq = 0.89, the symbol ★ to E/Eq = 
0.70, the symbol o to E/Eq = 0.77, the symbol o to E/Eq = 
0.80 and the line to E/Eq = 0.45.The filled symbols represent the 
second measurement realized for some E/Eq. 


the presence of flow, few large polyethylene spheres (Pe) 
were used as tracers in the surrounding fluids furnishing 
a qualitative visualization of the flow around the drop. 
Figure [5] illustrates the streamlines for the flow around a 
“clean” drop and the measured trajectories of two parti¬ 
cles around a drop covered with glass particles. For high 
surface coverage, the traces flow from the poles to the 
belt boundary. Tracers in the vicinity of the surface cov¬ 
ered by particles are still. Thus, the presence of particles 
suppress the electrohydrodynamic flow and limits it to 
the poles region (‘clean’ portion of the drop). 


B. Stronger DC electrical fields: drop tilt 

Particle-covered drops also exhibit Quincke-like (elec¬ 
trorotation) behavior. We find that the inclination angle 
does increase with field strength, but also exhibits un¬ 
steady behavior, see Figure[6] The appearance of nonzero 
tilt /3 (steady or unsteady) is defined as the threshold for 
rotation. 

Figure [6] (b) represents the time evolution of the tilt 
angle (3 (see Figure I] for the definition) measured for a 
drop covered of small glass spheres for different electrical 
fields. For low electrical fields, /? remains nearly zero, i.e., 
there is no tilt. The small angle is due to errors in the 
position of the electrical chamber relative to the camera. 
In stronger field the drop tilts. However, in contrast to 
‘clean’ drops which exhibit stationary inclination, the in¬ 
clination (/3) varies with time, see Figure [6j (b), similarly 
to the wobbling observed for drops with a low coverage of 
aluminum particles [24]. time scale of the transients is 
comparable to the duration of the observation, so we can 
not conclusively state that the behavior is periodic (e.g., 
oscillations). Figure [6] (a) represents the maximum in¬ 
clination fimax measured as a function of the electrical 
fields scaled by the critical value for Quincke rotation 
(E/Eq). Filled symbols correspond to repeated exper¬ 
iments. Below the rotation threshold, E c /Eq , Pmax is 
within the experimental error; E c denotes the critical 
electric field for the onset of electrorotation in our system 
(particle-coated drop). Above E c , Pmax increases with 
field strength E. Repeated measurements show that for 
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a given E , different inclinations are possible. Finally we 
observe sometimes that the a drop can spin while having 
no tilt. 

Figure [7} (a) shows the maximum inclination /3 max 
measured for drop covered with different particles as a 
function of E/Eq. The scatter of the data suggests that 
the tilting effect is strongly dependent on particles char¬ 
acteristics. For comparison, the solid black line repre¬ 
sents the theoretical dipole orientation prediction for a 
Quincke sphere Eq. [ 6 ] [33]. 



Pe 78%, 
a) 



G1 (7jic) 91% 



Ca= 1.27-E/E 0 = 0.83 Ca = 1.16 - E/E 0 = 0.894 Ca = 1.33 - E/E 0 = 0.855 


Cluster ejection Dynamic wings 


Implosion 



A1 88% A1 (3uc) 90% 


A1 (3 pc) 90% 




Ca= 1.71-E/E Q = 1.05 Ca= 1.71-E/E Q = 1.05 C&= 0.92 - E/E Q = 0.76 

Implosion Dynamic wings Dynamic wings 
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^3 d O 


FIG. 7: (a) Maximum inclination fimax measured over 4 minutes 
as a function of E/Eq. The vertical lines represent the onset of 
particle sustained motion for the different drops, and the black 
thick line represents the prediction for the dipole orientation for a 
solid particle given by Eq. [6] The symbols o, □, and v denote 
drops covered with A1 (r = 12 pm) with p = 93%, 88%, and 96% 
and a = 1.7mm, 1.8mm, and 2.2mm, the symbol <1 corresponds to a 
drop covered with Pe (p = 78%, a = 2.2mm), the symbols o, and a 
denote drops covered with G (r = 8.5 pm) with p = 100% and 93% 
and a = 1.2mm, and 2mm, the symbol + corresponds to a drop 
covered with G (r = 3.5/im) with <p = 100% and a = 1.7mm. For 
clarity, error bars are shown only for some drops as they are similar 
for all the measurements. (b) Electric field at which sustained 
drop tilt occurs as a function of the drop diameter. The filled 
symbol o represents a drop covered with Pe particles for p = 78%. 
The empty symbols represent drops covered with A1 particles, o 
for p = 95%, □ for c p = 93%, and A for p = 88%. The filled 
A symbol represents a drop covered with glass sphere (r = 3 pm) 
for p =. The filled □ and o symbol represent respectively drops 
covered with glass spheres (r = 12 pm) for p = 100% and p = 94%. 


Figure 7|(b) summarizes the threshold fields as a func¬ 
tion of the drop diameter for high surface coverage of 
different types of particles. Independently from the type 
of particle, the dimensionless critical field strength is in 
the range [0.5-0. 8 ]. i.a. two to four times lower than 
the thresholds measured for “clean” drops and around 
twice lower than for a solid sphere. Hence, high surface 
coverage of particles drastically reduces the critical fields 
strength for onset of rotation. 

For very strong electrical fields, the drops exhibit very 
peculiar behavior, see Figure [ 8 j (a), as drop implosion or 
bucking of the layer of particles leading to either ejection 
of cluster of particles similar to the ejection observed for 
lower surface coverage [24] , or the formation of ephemeral 
dynamic wings. Figure [ 8 j (b) illustrates ephemeral wings 
formation. 


FIG. 8: (a) Drop behaviors for very strong electrical fields, (b) 
Evolution of the “dynamic wings”. The scale bar corresponds to 1 
mm. 


IV. DISCUSSION 

In summary, our experiments show that surface- 
adsorbed particles suppress the electrohydrodynamic 
flow, enhance overall oblate drop deformation, and de¬ 
crease the threshold for rotational flow. Modification of 
the surface properties, e.g., surface tension, interfacial 
viscoelasticity and surface conductivity, could be respon¬ 
sible for these effects. Here we discuss each of these 
mechanisms and show that the oblate deformation of the 
particle-covered drop is well captured by a capsule model; 
the colloidal membrane endows the interface with shear 
elasticity which immobilizes the interface and suppresses 
the flow. The decrease in the tilt threshold, however, 
can not be explained at this time with the available shell 
model G3ST]. 


A. Modification of the surface tension 

The increase of deformation could be related to a lower 
effective interfacial tension, 7 e //, due to the particles. 
Fitting the data on Figure [4] with the Taylor model Eq. 
[5] yields 7 e ff in the range [0.12 - 0.16] mN/m, which cor¬ 
responds to a decrease of almost 97% relative to the ten¬ 
sion for the clean interface. Impact of surface absorbed 
particles on the surface tension has been studied in re¬ 
lation with Pickering emulsion and liquid marbles [42]. 
The presence of solid particles at the interface usually 
is associated reduction of the liquid/liquid interfacial 
tension mu mi. Theoretical modeling m shows that 
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for dense monolayer of mono-disperse spherical particles, 
7 e ff depends mainly on the contact angle of the particles, 
7e// = 7[1 — 7r(l — cosO) 2 / 4 a/3], and shouldn’t decrease 
below 50% of 7. This order of magnitude is confirmed by 
experimental measurement [43] [46] . For lower concentra¬ 
tion of particles, almost no variation of 7 were measured 
mm- Recent publications |HJ 48] on marble liquids 
point out the strong dependence of 7 e ff on size, concen¬ 
tration, surface free energy and water contact angle of 
the encapsulating powder. For non-densely packed par¬ 
ticles, the capillary interactions between particles can be 
neglected, and the “modified” surface tension can be ex- 
pressed as 7 eff = ^+{'y s i+lsa~l)A 0 /A gS], where 7 SL 
and 7 sa are respectively the surface tension at the pow- 
der/liquid and powder/air interface, and A and A$ are 
the total surface of the marble and the surface covered by 
particles. Estimate for all drops considered in our work 
show that 7 e ff should be relatively close to 7. Hence, 
a lower interfacial tension can not explain the strong in¬ 
crease of the drop deformation. Moreover, 7 e ff should 
vary with the electrical fields, as the packing changes with 
the electrical fields leading to a non linear variation with 
the deformation. In conclusion, a variation of the in¬ 
terfacial tension due to the presence of surface-adsorbed 
particles can not explain the strong deformation increase 
compared to a clean drop. 


B. Surface conductivity 

We can model the monolayer of surface-adsorbed par¬ 
ticles as a solid shell, with an equivalent surface conduc¬ 
tivity estimated as the particle conductivity divided by 
the particle diameter (which is a measure of the thickness 
particle monolayer), a s ~ cr p /(2r). The current conserva¬ 
tion boundary condition may be expressed as (assuming 
negligible charge convection, which is reasonable in weak 
fields and small drop deformations) cr ex E^ — (J- in E = 
—07 V s • E^. The solution shows that the deformation is 
described by the same discriminating function in the Tay¬ 
lor theory Eq. [5] but with R replaced by R + 2 R s where 
R s = (T s a/<T ex is the dimensionless surface conductivity. 
This shows that oblate deformation should be suppressed 
as the particle conductivity increases and eventually the 
drop deformation should change from oblate to prolate. 
In our experiments however we observe enhancement of 
the oblate deformation by the surface-adsorbed particles, 
not suppression. This implies that surface conductivity 
is unlikely to play a significant role in drop deformation. 


C. Interfacial viscosity 

Particles at interfaces behave as two-dimensional sus¬ 
pension. Accordingly, the increased dissipation arising 
from the particle motions results in increased surface vis¬ 
cosities f49H5T] . A drop with high surface viscosities ef¬ 
fectively acts as a drop with very high bulk viscosity. In 



FIG. 9: Drop deformation as a function of the capillary 
number. Symbols the same as in Figure [4] The solid line 
represents the Taylor theory Eq. [5] for A = 0.07 , the 
dashed line is Eq. [8] for A —>> 00, and the lines for capsule 
model Eq. [9] with shear elastic modulus 1.437. 


the limit of high viscosity ratio Taylor’s law gives 


D = 


9 (-19+ (5 + 9R + 5R 2 ) S) 
° a 805(2+ R) 2 1 
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Increased surface viscosity also suppresses the surface 
flow, see Eq. [4j 

Figure [9] shows the comparison between the experimen¬ 
tal data and the theoretical curve for an infinitely viscous 
drop Eq. [8j The theoretical curve is close to the data spe¬ 
cially for low capillary numbers but still underestimates 
the measured deformation. Hence, the effect of the pres¬ 
ence of particles can not be completely explained by the 
suppression of the electrohydrodynamic flow. Moreover, 
we also observe that drop deformation continues to in¬ 
crease even above the rotation threshold (see Figure [4]) 
in contrast to the theoretical prediction for high-viscosity 
drops [40] . 


D. Interfacial elasticity 

A monolayer of densely packed particles at a liquid in¬ 
terface behaves similarly to an elastic sheet [521456] and 
can support anisotropic stresses and strains: it buck¬ 
les in uniaxial compression and cracks under tension. 
This solid-like behavior can be characterized in terms of 
Young’s modulus and Poisson ratio derived from simple 
theoretical considerations see [S2 El EH- Recent the¬ 
oretical model [58] shows that presence of a significant 
interfacial shear elasticity (greater than the interfacial 
tension) and of a nonzero interfacial compression elastic¬ 
ity can cause drops to buckle. In our experiments, the 
particle layer also eventually exhibit buckling, see Figure 
[SJsuggesting that the particle-coated drop behaves like a 
capsule (a drop enclosed by an elastic membrane). 
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The deformation of capsules in electric fields has been 
considered in only two studies mm- Ha and Yang m 
developed a small deformation theory similar to the one 
for a capsule in external linear flows [6T . Since the par¬ 
ticles at the interface are closely packed, the monolayer 
of jammed particles can be much more easily sheared 
than compressed (compacted). Accordingly, the parti¬ 
cle monolayer can be modeled as an incompressible two- 
dimensional material and drop deformation is given by 
(see Appendix for details) 

= C«1 
(9) 

where G is the elastic shear modulus. Figure [9] show the 
data fit with y/G = 0.7. The value of G is in reasonable 
agreement with value deduced from the Young’s modulus 
of particle rafts [52] . 
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E. Threshold for drop tilt and electrorotation 


stronger fields, drops tilt due to the Quincke effect but 
the onset is significantly lowered compared to the clean 
drop. The decrease can not be explained with the ex¬ 
isting spherical shell model and likely due to drop as- 
phericity. Even stronger electrical fields give rise to more 
exotic behaviors as ephemeral dynamic wings or drops 
implosion. The similarity between the particle-covered 
drop and a capsule provides promising new insight into 
the impact of particles on interfacial dynamics. Our 
findings open questions ranging from understanding the 
fluid/solid transition that occur for high surface coverage 
to stability of Pickering emulsions. We hope our work will 
stimulate further research on the electrohydrodynamics 
of particles at interfaces. 
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To estimate the effects of the particle monolayer on the 
onset of electrorotation, we consider the spherical shell 
[371 ED- We can model the particle-covered drop as a 
layered sphere with effective properties e v and a v . For a 
thin shell [62] 


^ in 
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where the shell thickness corresponds to the particle di¬ 
ameter 2r. While S is less affected by the particles, R 
can be significantly changed due to the wide range of 
particle conductivities. For example, Eq. [To] shows that 
presence of a resistive (very low-conductivity) layer de¬ 
creases R and thereby lowering the rotation threshold 
Eq. [7] However, since this decrease is also observed for 
conductive particles, surface conduction most likely not 
a relevant explanation. A possible explanation may be 
drop asphericity; prolate rigid ellipsoids do have lower 
threshold for electrorotation than spheres [63] [64 ] . This 
issue requires further examination both experimentally 
and theoretically. 


V. CONCLUSIONS 


Appendix A: Small-deformation theory for an 
initially spherical elastic capsule placed in a uniform 
DC electric field 


An uniform DC electric field, E = Eoz, exerts elec¬ 
tric pressure and shear on a sphere with non-capacitive 
interface (continuous electric potential) [5Tj 


t ei = p e[ (1 + 3 cos 20) r + rj 1 sin 2 00 (Al) 


where 


P el = _ 3 (-ft 2 + 1 — 2/5) 


4(ft + 2) 2 


HR-m 
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The electric tractions deforms the interface; the material 
particles at the interface move to new position described 
by 


x s 


= a 


(1 + s (1 + 3 cos 20)) r + us\n20Q 


(A3) 


The in-plane displacement does not generate overall 
shape change but creates elastic stresses that oppose the 
electric shear and immobilize the interface. The devia¬ 
tion from sphericity is qualified by s 


We experimentally study the effect of high concentra¬ 
tion of surface-adsorbed particles on drop deformation 
in a uniform DC electrical field. The fluid system con¬ 
sists of a silicon oil drop suspended in castor oil, both 
very weakly conducting liquids. A broad range of par¬ 
ticle sizes, conductivities, and shapes is explored. In 
weak electric fields, the presence of particles enhances 
the deformation compared to a particle-free (clean) drop 
and suppresses the electrohydrodynamic flow. Drop de¬ 
formation is well described by a capsule model, which 
treats the particle monolayer as an elastic sheet. In 


D = 3s , d|| = r s ( 0) = 1 + 4s , d± = r s (ic/2) = 1 — 2 s 

(A4) 

For small deformations the membrane behaves as an 
elastic-Hookean material and the elastic stresses are 

[S3 E3 [65] 

t m = (1 + 3 cos 20) f + r s m sin 200 , (A5) 

where 

r m = 2 (2 Ca~ 1 u + 3Ca~ 1 (2s + u)) p m = 2Ga“ 1 (2s+'u). 

(A6) 












Ca s and Ca a are the dimensionless shear and extensional 
elasticities; Ca s = e^E^a/G. Balancing the electric and 
elastic stresses, Eq. |Al| and Eq. |A5[ yields 

s = Ca s - (3 Ca s +2Ca a )^~ (A7) 

o o 


In the limit of inextensible membrane, Ca a 0, we ob¬ 
tain Eq. [9] 
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